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ABSTRACT 

The  phenomenon  of  thermal  lensing  was  investigated  in  water  using  a  Z-scan  method  and  corresponding  first-order 
mathemancal  models.  Data  from  first-order  thermal  lensing  models  and  AJBCD  beam  propagation  methods  were  used  to 
simulate  the  non-lmear  absorption  of  water  held  in  a  thin  sample  cuvette  for  a  Z-scan  optical  set  up  of  CW  cases  at  1313 
nm.  The  single  beam  closed  aperture  Z-scan  was  then  used  to  determine  the  non-linear  absorption  at  1313  nm  for  water  in 
10  mm  and  2  mm  cuvettes  at  48.00,  16.80,  9.80  and  2.83  mW  then  compared  to  the  first-order  model  data.  The  results 
from  the  closed  aperture  Z-scan  were  also  used  to  back  calculate  the  spot  size  in  the  far  field  for  comparison  to  the  model’s 
prediction  of  the  beam’s  temporal  response.  Experimental  Z-scan  data  were  found  not  to  correlate  strongly  with  our  first- 
order  model  suggesting  the  need  for  higher  order  models  to  successfully  predict  spot  size  in  absorbing  media  inside  the 
Rayleigh  range. 
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1.  INTRODUCTION 

Damage  thresholds  are  important  for  providing  safety  guidelines  in  the  American  National  Standards  Institute,  ANSI.  As 
a  result,  these  safety  guidelines  are  used  in  areas  which  employ  the  use  of  laser  technology.  For  example,  retinal  laser 
surgeries  are  used  to  treat  ocular  disorders  such  as  macular  edema  and  central  serous  retinopathy  by  coagulating  small 
areas  of  abnormal  micro  vascular  activity/  s  Predicting  safe  exposure  levels  and  subsequent  damage  thresholds  for  pulsed 
or  continuous  wave,  CW,  lasers  for  treatment  wavelengths  at  a  given  power  require  knowledge  about  the  spot  size  of  the 
beam  on  the  retinal  plane.  If  beam  parameters  such  as  spot  size  and  power  at  the  retina  are  known  then  a  damage  threshold 
can  be  derived  based  on  tissue  properties  and  the  fluence  rate  on  the  retina.4, 5  It  is  hypothesized  that  thermal  lensing 
within  the  human  eye  yields  a  significant  amount  of  change  in  beam  geometry  impacting  the  resulting  damage  threshold 
value  for  a  variety  of  spot  size  conditions. 

Thermal  lensing  is  a  well  known  phenomenon  in  analytical  chemistry  where  a  temperature  gradient,  AT,  forms  in  an 
absorbing  media  resulting  in  a  change  in  density,  p,  of  the  material  where  absorption  occurred  leading  to  the  change  in 
index  of  refraction,  n,  with  respect  to  temperature.6  *  This  effect  yields  a  thermo-optic  coefficient  represented 
mathematically  as  dn/dT.  Whenever  the  index  of  refraction  changes,  the  path  the  beam  takes  through  the  media  is  altered 
from  the  original  path  resulting  in  a  shift  m  focus.  In  media  which  exhibit  a  decrease  in  density  as  temperature  increases, 
such  as  ocular  tissue,  the  index  of  refraction  decreases  with  respect  to  the  temperature  rise,  a  negative  dn/dT,  resulting  in  a 
negative  lens  like  effect  also  known  as  thermal  blooming.  A  negative  lens  effect  would  also  mean  the  focal  point  w'ould 
shift  behind  the  original  focal  plane. 

The  long  term  research  goals  in  this  area  are  to  develop  thermal  lensing  models  for  ocular  media  derived  from  earlier 
thermo-optic  models  developed  for  solid  state  optics  research.  9’ 1C  11  These  adapted  thermal  lensing  models  will  be 
compared  to  experimental  data  for  both  thick  and  thin  samples  for  validation,  then  used  to  provide  information  to  retinal 
damage  models  created  by  Thomas  et  al.  m  2003. 4’ 12  Ultimately,  this  research  will  impact  the  way  in  which  predictions 
and  assessments  for  damage  on  the  retina  are  made  in  order  to  provide  critical  information  to  the  medical  community  for 
eye  surgeries  and  retinal  imaging  techniques  which  utilize  laser  technologies.  In  this  study,  a  first-order  thermal  lensing 
model  was  created  in  Lab  VIEW™  (National  Instruments)  to  simulate  the  thermal  lensing  properties  of  the  human  eye. 
First-order  model  data  were  compared  to  experimental  data  collected  using  a  single-beam  closed  aperture  Z-scan  set  up.6 
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2.  METHODOLOGY 


2.1  First-order  model 

In  the  1960’s,  Gordon  et  ah*  discovered  as  media  absorbs  energy,  the  pathlength  of  light  through  a  media  changes, 
resulting  m  changes  of  a  beam  profile  within  a  sample.  Gordon  et  al.’s  *  discovery  became  known  as  thermal  blooming  or 
thermal  lensmg  since  many  media  exhibited  characteristics  of  broadening  the  beam  upon  exposure  to  light  energy.  To 
successfully  model  thermal  lensing,  time  dependent  terms  for  describing  both  thermal  effects  and  beam  geometry  within  a 
sample  must  be  determined. v  1  1 

As  a  media  absorbs  energy,  physical  characteristics  such  as  density  and  temperature  change.  These  physical  changes 
impact  the  index  of  refraction,  n.  Van  Stryland  et  al.11  describe  n  m  terms  of  the  change  m  density,  p,  and  temperature  T 
by  a  term  for  the  thermal  expansion,  /  dp  and  a  term  for  other  temperature  changes,  dn/dT  .  The  term  for  thermal 
expansion  is  commonly  referred  to  as  electrostriction. 
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The  electrostriction  component  can  be  approximated  by  a  constant  if  the  time  it  takes  for  an  acoustic  wave  to  travel  the 
distance  of  the  minimum  beam  spot  is  less  than  the  time  of  the  pulse  w  idth  or  pulse  on  time,  i.11  This  would  hold  true  for 
cases  involving  long  pulse  durations  or  contmuous  wave,  CW.  This  simplifies  Eq.  1  to  the  thermo-optic  effect  m  Eq.  2. 

To  consider  the  thermal  effects  for  the  first-order  model  a  Gaussian  beam  in  TEMqo  mode  w  as  assumed  to  propagate 
through  an  absorbing  media  where  the  temperature  along  the  beam  path  would  be  greatest  in  the  center  and  coolest  at  the 
outer  edge  of  the  beam  profile.  This  would  mean  the  index  of  refraction  in  areas  at  higher  temperatures  would  be  lower 
than  the  index  of  refraction  in  areas  of  cooler  temperatures,  ultimately  describing  how  a  media  can  exhibit  a  negative  lens 
like  behavior.  The  thermal  profile  within  an  absorbing  media  was  described  by  the  radius,  r  ,  from  the  optical  axis, 
position,  z,  along  the  optical  axis,  and  time,  t,  resulting  in  Eq.  2.9, 10 
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Where  n©  is  the  background  index  of  refraction.  Electrostriction  was  neglected  for  the  purposes  of  this  study  as  pulse 
widths  were  greater  than  the  acoustic  transit  times  for  all  cases  represented  here. 

The  mathematical  model  used  for  the  thermal  lensing  effect  originates  from  the  Green’s  Function  solution  to  the  thermal 
diffusion  equation  (Eq.  3).  13 

V2r(F,0--arVV)  =-4*S(r,0  (3) 

7]  dt 

Where  T  is  temperature,  t  is  time,  r  is  the  radial  vector,  rj  is  the  thermal  diffusivity,  and  S  is  the  source  term. 

Assummg  a  cylindrical  symmetric  medium  where  diffusion  along  the  optical  axis  is  zero  (dz/c)T  =  0),  then  the  solution 
to  the  diffusion  equation  at  radius,  F  and  time,  t,  can  be  described  by  9 
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Where  G{r  ,t\  r'  ,t')  =  0  for  all  positions  where  r  is  on  the  surface  of  the  cylindrical  media  and  G(r,f;  r\t')  =  0  for  all 
times,  t  <  t\ 

The  temperature  rise  in  the  medium  at  position  r  can  be  determined  by  integrating  over  all  r '  in  the  medium  In  order 
to  accurately  represent  the  accumulating  thermal  effects  of  a  continuously  pulsed  or  chopped  source,  the  temperature  rise 
in  the  media  will  be  related  to  the  current  pulse  and  any  residual  contributions  from  previous  pulses.1  ‘ 

dt'  ^  dr'S(r,)G(f,r,r,n  +  Yj  j~  dr'S(r')G(r,t;r\t'n)  (5) 

n=l  " 


Where  tm  is  the  time  during  the  mLi  on  pulse.  An  =  (tm  +  x’)+(n  -  1)(T  +  t’)  with  T  =  on  time  and  x’  =  off  time. 
The  temperature  change  for  when  the  pulse  is  off  can  then  be  described  by  the  events  following  the  m^  pulse  b 

T(r,t'm  ;m)  =  £  jf dt' dr' S(f')G(r,t;r' ,t\ n )  (6) 

n- 1  m/1 

Where  Cn  =  t*m  +  (n  -  1)(t  +  r).  Fig.  1  depicts  the  timing  diagram. 
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Figure  1  The  timing  diagram  used  in  the  thermal  lensing  model  for  m  pulses  with  on  time,  x,  and  off  time,  x\  (a)  and 
(b)  depict  the  sequences  used  for  when  the  observation  point  is  taken  when  the  pulse  is  on  and  off  respectively.  Note 
the  point  of  observation  is  actually  the  zero  point  in  time.  Equations  for  An  and  Q  are  described  in  Eq.  5  and  6  t^  is 
the  time  from  the  last  rising  edge  of  the  most  recent  on  pulse  to  the  point  of  observation  and  t^’  is  the  time  from  the  last 
falling  edge  of  the  most  recent  on  pulse  to  the  pomt  of  observation  taken  while  the  pulse  is  off.  (Figure  adapted  from 
Swofford  et  al.)10 

Using  the  formulism  by  Fang  and  Swofford,  Eq.  2  can  be  re-written  as  10 
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To  describe  the  temporal  dependency  for  the  thermal  response,  an  expression  for  the  source  term  can  be  found  by  applying 
Beer’s  law  for  the  heat  generated  from  absorption  per  unit  length  through  the  sample  and  assuming  a  Gaussian  beam  '  1 
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Where  p*  is  the  linear  absorption  coefficient,  0)  is  the  1/e2  beam  radius,  k  is  the  thermal  conductivity  and  P0(t)  is  the  power 
as  a  function  of  time,  t.  Cases  using  a  square  pulse  can  allow  P0(t)  to  be  a  constant. 

Cases  involving  one  or  more  square  temporal  pulses  can  be  written  in  terms  of  the  on  and  off  portions  of  the  pulse  shown 
in  Eq.  9  and  10  respectively. 
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Where  is  the  linear  absorption  coefficient,  k  is  the  thermal  conductivity,  P  is  the  power  when  the  pulse  is  on,  0)  is  the 
1/e2  beam  radius  (assumed  to  be  constant  over  a  very  small  increment,  A Z  ),  tm  is  the  time  from  the  last  rising  edge  of  the 
most  recent  on  pulse  to  the  point  of  observation,  X  is  the  pulse  on  time,  and  k  is  the  critical  time  defined  by  tc  =  T|/8cOo. 

To  propagate  the  Gaussian  beam  through  space,  we  used  the  AJBCD  law.7  14  The  Gaussian  beam  is  described  by  a 
complex  number  at  spatial  position,  z,  w  ith  the  radius  of  curvature,  R,  representing  the  real  portion  and  the  beam  radius 
and  wavelength  for  the  imaginary  component. 
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Yariv  details  how  to  model  a  Gaussian  beam  moving  through  a  diverging  lens-like  media  described  by  a  quadratic  index 
of  refraction  by  the  following  transformation  matrix. 
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Where  Az  is  a  small  spatial  increment  along  the  optical  axis  and  X  is  defined  in  Eq.  13. 10 
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The  new  complex  beam  parameter  is  updated  by  the  previous  value  of  q(z)  and  the  matrix  in  Eq.  14  by  the  following 
relationship. 


Aq(z)M+B 
Cq(z)old  +  D 


(14) 


The  solution  is  found  by  iteratively  calculating  the  beam  radius  as  a  function  of  spatial  position,  0)(z),  through  the  sample 
using  an  initial  guess  for  the  case  of  no  thermal  lensing.  This  process  yields  a  self-consistent  description  of  co(z)  at  any 
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observation  time,  t.  The  resulting  data  for  co(z)  can  then  be  used  to  evaluate  the  thermodynamic  interaction  of  the  beam 
with  the  sample. 

2.2  Single  beam,  closed  aperture  Z-scan  experiment 

Experimentally  characterizing  thermal  lensing  in  a  media  can  be  accomplished  by  using  the  Z-scan  method.6  1S' 16, 1  In  a 
Z-scan,  the  relative  transmittance  is  read  by  a  photodiode  placed  behind  an  aperture  as  the  sample  is  moved  in  incremental 
steps  in  front,  through  and  behind  the  depth  of  focus  of  a  focusing  lens.  Optical  elements  such  as  the  focusmg  lens  and 
neutral  density  (ND)  filers  are  placed  before  the  chopper  in  order  for  any  thermal  lensing  effects  within  optical  elements  to 
reach  steady  state.  Resulung  data  are  measurements  of  relative  beam  intensity,  related  by  transmittance  detected  over  the 
surface  area  of  the  photodiode,  where  a  value  for  dn/dT  can  be  found  by  determining  the  best  fit  to  the  experimental  Z-scan 
data. 
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Figure  2:  The  Z-scan  set-up  used  for  the  experimental  data  collection.  Sample  cuvettes  were  2  and  10  mm  thick  and 
contained  DI  FLO.  Positions  to  the  left  (in  front)  of  the  focal  point  of  the  lens  are  negative  and  positions  to  the  right 
(behind)  of  the  focal  point  are  positive.  All  transmittance  values  were  normalized  to  the  transmittance  detected  when 
the  sample  was  placed  in  the  beam  path  outside  the  depth  of  focus  of  the  lens 

In  a  closed  aperture  z-scan,  the  photodiode  is  placed  in  the  far  field,  in  this  case  it  wras  placed  117  cm  from  the  pump  beam 
focusing  lens.  The  probe  beam  transmittance  was  then  normalized  relative  to  the  energy  detected  when  the  shutter  was 
open  and  the  sample  was  70  cm  after  the  focus;  a  position  where  no  thermal  lensing  should  occur.  Upon  exposure  to  the 
beam,  the  thermal  lensing  effect  for  a  negative  lens  like  media  will  result  in  an  increase  and  decrease  in  relative 
transmittance  for  cases  mside  the  Rayleigh  range  before  and  after  the  focus  respectively.  The  sample  was  moved  along  the 
z-axis  to  achieve  transmittance  readings  for  various  areas  within  the  Rayleigh  range.  The  inflection  point  is  observed  at  the 
focus.  An  example  of  what  closed  aperture  z-scan  data  appear  as  for  a  negative  lens  case  is  shown  in  Fig.  3  below. 
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Figure  3:  An  example  of  typical  data  from  a  closed  aperture  z-scan.  Z-position  refers  to  the  position  of  the  sample 
relative  to  the  focus  marked  as  zero. 

A  Crystal  Laser  Model  IRCL-150-1313-P-L  was  used.  A  Galilean  telescope  was  configured  to  create  a  collimated 
expanded  beam  with  1/e2  diameters  of  3.5  mm.  A  focusmg  lens  of  400  mm  made  of  BK7  glass  was  used.  A  knife  edge 
test  was  conducted  to  determine  the  focal  position  of  the  lens.  The  pump  beam  focusing  lens  was  placed  before  the 
chopper  to  allow  for  any  thermal  effects  within  the  BK7  glass  to  achieve  steady  state  throughout  the  duration  of  the 
experiment. 
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The  sample  cuvette  containing  de-ionized  ultra-filtered  water  was  moved  along  the  z-axis  by  a  large  micrometer  mount. 
Blank  measurements  for  an  empty  cuvette  were  determined  to  yield  no  change  in  normalized  transmittance,  thus  assuring 
thermal  lensing  of  the  cuvette  for  the  pump  wavelength  was  negligible  in  this  experiment. 

The  power  of  the  beam  was  lowered  by  using  a  Thorlab  NEK01  neutral  density  filter  set.  When  a  ND  filter  of  0.5,  1 .0,  2.0 
and  3.0  were  inserted  into  the  pump  path  a  power  of  48,  16.8,  9.8  and  2.83  mW  were  delivered  to  the  sample  respectively. 
The  system  was  re-aligned  each  time  a  new  ND  filter  was  added  or  removed  from  the  optical  set  up.  The  photodiode  was  a 
United  Detector  Technology  8028  and  was  connected  to  a  Melles  Griot  large  dynamic  range  transimpedance  amplifier 
where  the  signal  was  then  read  into  a  Labview  program.  Data  were  then  processed  in  Matlab  to  compare  to  the  first  order 
model. 

The  aperture  size  was  selected  to  achieve  readings  between  30-75%  of  the  maximum  power  read  at  the  selected  detector 
position.  The  aperture  diameter  was  2.5  mm  for  each  experiment.  The  z-positions  ran  from  -1 8  to  52  cm  in  incremental 
steps  of  1  cm.  The  m  value  for  the  beam  profile  was  found  to  be  1.16. 

Experimental  cases  were  investigated  with  the  following  z-scan  parameters  summarized  in  the  Tables  below. 
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wavelength 

(nm) 

Time 

between 
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(s) 

Sampling 
rate  (Hz) 

Measured  lens 
focal  length  (using 
knife  edge)  (cm) 

1/e2  Beam  radius 
(at  lens)  (mm) 

Distance  from 
lens  to  detector 
(cm) 

Iris 

diameter 

(mm) 

1313 

20 

2500 

42 

1.74 

117 

2.5 

(b) 


Exposure 
duration  (sec) 

Sample  thickness 
(mm) 

Pump  power  at 
sample  (mW) 

Temperature  of  room  (sample 
at  room  temp)  (°C) 

0.1 

2 

2.83 

18.3 

9.8 

17.3 

16.8 

16.1 

48 

16.1 

10 

2.83 

18.3 

9.8 

17.3 

16.8 

16  1 

48 

16.1 

(c) 


Sample  thickness 
(mm) 

Exposure 
duration  (sec) 

Pump  power  at 
sample  (mW) 

Temperature  of  room  (sample 
at  room  temp)  (°C) 

10 

0.5 

2.83 

18.9 

48 

206 

1 

2.83 

18.9 

48 

20 

Table  1 :  Single  beam  closed  aperture  z-scan  experimental  parameters  in  this  comparison  study,  (a)  Standard  settings 
used  in  each  experiment  conducted,  (b)  Parameters  used  for  comparison  of  varying  power  and  sample  thickness,  (c) 
Parameters  used  for  comparison  of  varying  exposure  duration. 

3.  DATA  &  RESULTS 

Several  studies  have  been  conducted  on  the  linear  absorption  coefficient  for  water.1*’ 19  For  the  purposes  of  this  study,  the 
absorption  coefficient  data  from  Hale  et  al.18  were  used  to  approximate  an  a  for  1313  nm  yielding  126  1/m.  The  value 
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used  for  k  was  0.6  J/(msk).  Determining  a  dn/dT  for  13 13  nm  was  done  by  fitting  a  polynomial  to  data  published  by 
Schiebener  et  al.  22  The  background  index  of  refraction  for  water  at  1313  nm  was  1.315  and  was  assumed  to  be  for  an 
ambient  temperature  of  25  °C. 

The  first  order  model  was  run  for  each  of  the  expenmen tal  equivalents  described  in  Table  1 .  The  results  are  shown  in 
Fig’s.  4  and  5  below. 


(a) 


Singi«  B**m  Z-Sc»r\.  100m*  Smjl*  Puk**.  Pow»r  K  2  83mW 


(b) 


Smgw  B«*m  ZSc*n.  100m*  Pul**  Po**»r  ■  9  8mW 


(C) 


Smgl»  Beam  Z  Scan  IQOma  Sngla  Pvlaa,  Powa»  ■  16  8mW 


(d) 


Smgia  B*am  2- Scar  IQOma  Smgt*  Pda*  Powaf  »  48mW 


Figure  4-  Results  of  the  model  and  expenment  for  2  mm  (circles)  and  10  mm  (x’s)  with  a  100  ms  exposure  for  (a)  2.83 
mW,  (b)  9.8  mW,  (c)  16  8  mW  and  (d)  48  mW. 

Modeling  results  deviated  further  from  the  data  as  sample  thickness  increased  and  power  increased.  The  data  for  the  2  mm 
thick  sample  in  Fig.  1(a)  was  very  noisy  and  difficult  to  detect  even  with  several  scans  and  multiple  averaging  techniques. 
The  model  seems  to  fit  the  2  mm  sample  more  successfully  in  three  of  the  four  power  settings  compared  to  the  10  mm 
sample. 


Smgi«  B«*m  Z-Sc*n,  10mm  S*mpl«.  Pow»r*  iflmW 


Singi*  B«*m  Z  Sc*n.  10mm  S»mpi* .  Fow*r  *■  2  B3mW 


Figure  5:  Results  of  the  model  and  expenment  for  10  mm  at  0. 1,  0.5  and  1 .0  s  exposures  for  (a)  48  mW  and  (b)  2  83 
mW. 


The  modeling  results  shown  in  Fig.  5  reflect  how  the  model  deviates  from  the  data  further  as  the  length  of  exposure 
increases.  Fig.  5  also  illustrates  how  the  model  over  estimates  the  thermal  lensing  effect  further  as  power  increases. 
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4.  DISCUSSION 


Originally,  the  model  was  used  to  try  to  find  a  dn/dT  which  fit  the  data.  While  running  simulations  with  the  model,  it  was 
discovered  small  increases  in  the  magnitude  of  the  linear  absorption  coefficient,  a,  resulted  in  an  increase  in  the  maximum 
inflection  point  before  the  focus,  causing  a  greater  exaggeration  m  the  prediction  of  the  thermal  lens  effect.  Several  values 
of  a  were  investigated  from  a  variety  of  literature  sources  including  Hale  et  al  and  Kou  et  al.18, 19, 2"  The  smallest  value  for 
a  seemed  to  fit  the  single  beam  z-scan  more  closely,  however,  a  was  not  the  only  variable  which  yielded  differences  in  the 
modeling  results. 

Studies  investigating  the  index  of  refraction  at  specific  temperatures  for  water  at  1.313  Jim  could  not  be  found,  though 
Harvey  et  al.21  suggests  interpolating  n(T)  from  Schiebener  et  al.20  is  a  good  approximation  at  1 .313  Jim.  The  interpolation 
procedure  used  to  approximate  n(T)  is  heavily  dependent  on  a  known  value  of  n  at  one  known  temperature  to  be  used  as 
the  reference  n  from  which  the  approximation  is  made.  The  approximation  was  only  done  for  data  between  10  and  90  °C 
at  0.1  MPa.  This  background  index  of  refraction  influenced  the  intensity  of  the  thermal  lensing  effect.  If  the  value  for 
dn/dT  was  held  constant,  then  the  modeling  result  yielded  little  to  no  change  in  the  relative  intensity.  Changes  in  the 
polynomial  coefficients  used  to  describe  the  function  n(T)  in  our  model  as  determined  from  Schiebener  et  al.20  data 
resulted  in  over  or  under  predicting  the  thermal  lensing  effect  as  the  slope  of  dn/dT  increased  and  decreased  respectively. 

Changing  the  value  of  k  influenced  the  modeling  results  as  well.  Increasing  k  by  as  little  as  1  %  yielded  a  noticeable 
decrease  in  the  peak  inflection  point  before  the  focus.  The  value  of  K  used  m  this  study  was  taken  from  the  CRC 
Handbook  for  water  at  298  K.  2  Future  thermal  lensing  models  may  need  to  consider  k(T). 

As  a  result  of  the  sensitivity  of  the  first  order  model  to  a,  k  and  dn/dT,  the  conclusion  was  made  to  not  use  the  model  to 
find  dn/dT  by  fitting  to  the  data.  Since  the  modeling  results  indicated  the  trends  of  over  exaggerating  the  thermal  lensing 
effect  for  thicker  samples  as  well  as  for  longer  exposure  durations,  the  selection  of  variables  like  K  and  a  could  become 
more  important  for  the  long  term  goal  of  developing  a  model  to  handle  thick  samples  like  where  strong  thermal  lensing 
effects  occur  such  as  in  the  human  eye.  Instead  further  investigation  into  the  selection  and  sensitivity  of  the  model  to  a 
and  k  will  be  done  for  future  models. 

Given  the  trends  seen  in  Fig’s.  4  and  5  of  the  model  fitting  lower  power  and  thinner  sample  cases  more  closely,  one  would 
predict  the  case  for  the  2  mm  thick  sample  at  2.83  mW  would  give  the  best  fit.  Unfortunately,  the  results  from  the  single 
beam  z-scan  data  at  lower  powers  (Fig.  4a)  were  inherently  noisy  making  the  determination  of  the  quality  of  the  fit 
difficult.  Further  data  collection  for  lower  powers  needs  to  be  conducted  in  future  studies. 

Preliminary  experimental  data  indicate  the  first-order  model  is  not  an  appropriate  approximation  for  thermal  lensing  for 
thick  sample  cases  with  strong  thermal  lensing  effects.  It  is  likely  the  aberrations  generated  in  the  thermal  lens  experiment 
actually  break  the  assumption  of  a  parabolic  distribution  for  n(r)  as  described  in  Eq.  1 1 .  Higher  order  models  are  under 
current  investigation.  However,  the  experiment  needs  to  be  repeated  under  reduced  noise  levels  within  the  instrumentation 
to  see  if  sensitivity  to  lower  power  levels  can  be  achieved  to  more  accurately  test  the  limitations  of  the  first  order  model. 
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